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BOUNDARY VALUE AND EXPANSION PROBLEMS; FORMULA- 
TION OF VARIOUS TRANSCENDENTAL PROBLEMS.* 

By R. D. Caemichael. 

Proceeding under the guidance afforded by a previous investigation f 
of the algebraic basis of the theory of certain boundary value and expansion 
problems we here develop the formal aspects of several transcendental 
expansion problems which are rather different externally but are abstractly 
of much the same essence. In each case the work is carried forward to the 
point of derivation of the form of the expansion and the values of the 
coefficients; and the problem of determining the nature of the functions 
having expansions of each of the kinds enumerated is thus clearly set forth. 
It seemed desirable to present in a single paper this formulation of several 
problems, so that the reader may most conveniently compare their elements 
of similarity or difference; and to leave to further investigation the detailed 
solution of each problem deserving such further treatment. It will be 
apparent to every one that we have given merely a selection of problems 
from a much larger number which emerge in a similar way from the con- 
sideration of various sorts of limiting cases of the algebraic systems dealt 
with in the paper already mentioned. 

In § 1 we treat the boundary value and expansion problem arising in 
connection with a single linear differential equation involving a single 
parameter. In § 2 the corresponding problem for a system of first order 
equations is more briefly treated and the results are carried over by 
Volterra's limiting process so as to yield the formulation of a similar expan- 
sion problem for certain integro-differential systems; and in § 3 the results 
of § 2 are rapidly generalized to problems involving v parameters. Various 
problems involving difference and integro-difference equations are treated 
in § 4; these are analogous to the corresponding ones for differential and 
integro-differential equations developed in §§ 1, 2, 3. In § 5 is given a 
brief statement of a great variety of mixed problems involving one or more 
parameters. The character of the problems arising from partial differential 
equations and integro-differential equations is developed rapidly in § 6 
and certain typical cases are briefly treated separately. A particular case 
of problems arising from partial differential equations is one which is con- 
nected with the theory of vibrating plates; this is briefly treated in § 7. 

* Presented to the American Mathematical Society, April, 1920. 
t American Journal op Mathematics, Vol. XLIII (1921), p. 69. 
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Cabmichael: Boundary Value and Expansion Problems. 233 

More details are given in § 1 than elsewhere in the paper, the justification 
being that those which are omitted are fairly obvious to one who has in 
mind the corresponding developments in § 1. 

1. The Problem Arising from a Single Differential Equation with One 
Parameter. — In connection with the linear differential expressions 

m s ln d^ + l ^^ + • • * + h Tx + loU ' 

T , , d m u , d m-1 w . , du . 

l(M) ~ 9m dx^ + 9 ^ 1 d^ i+ " ' + 9l 'dlc + 9oU ' m < n ' 

let us consider also the Lagrange adjoint expressions 

d n v . . dv , 

= m n - r -+ ■■■ + mi — + m v, 
dx n dx 

Mi(v) ^ hm d^ + '" + hi d^ +hoV - 

Here it is understood that m may- be zero so that Li(u) may have the par- 
ticular value gou. We assume that the coefficients h, gk, fn-k, hk, with 
their derivatives of all orders, are real-valued single-valued and continuous 
functions of x in an interval (ab) and that l n and g m do not vanish in (ab) . 
Without loss of generality we take l„ to be positive in (ab). With the linear 
differential equation 

(1) L(u) + \Liiu) = 

and certain n linear homogeneous conditions in u(a), u'(a), • • • , w (n-1) (a), 
u(b),u'(b), •••,«<«-»&, 

(2) tfi(w) = 0, U 2 (u) = 0, •••, U n (u) = 0, 
we associate the adjoint linear differential equation 

(3). M(v) + XJfi(») = 

and n like adjoint conditions 

(4) V n+1 (v) = 0, V n+2 (v) = 0, • • •, V 2n (v) = 

in such wise as to generalize the corresponding problem treated by Birkhoff * 
for the case when Li(u) = g u. 

We must first make precise the boundary conditions. In doing this 
as well as in developing the other parts of this section we follow the sugges- 

* Transactions of the American Mathematical Society, 9 (1908), 373-395; 219-231. 
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tions afforded by Birkhoff's treatment of the special case. If we write 

R(U, V) = M£5+ [^ ~ ^ (LV) ] £5 
+ 

+ [*-&™+£*<m -■■■+<- ^^H* 

fli(«, »)* = [^]^^+ • • • + ygit - ^ (<^) + ^ (go) 

+ (-D"- 1 ^i(^)]^ 

then it is well-known and easily verified that we have the identities 

(5) vL(u) — uM(v) = -7- {fl(w, v) }, vLiiu) — uMi{v) = -=- {Ri(u, v) } . 

Let us now consider the quantity 

[fl(«, v) + Xldiu, v)J= h a 

as a bilinear form in the two sets of variables w (n-1) (a), w (n_2) (a), • • ■ ,u'{a), 
«(o), «<"-»(6), m ( "~ 2) (6), • • •, u'(6), «(&) and »(o), »'(o), • • •, « ( "-«(o), »(6), 
«'(&), • • •> « (n-1) (&), and let the form be arranged in a square array whose 
columns contain the u's in the order written and whose rows contain the 
v's in the order written. The matrix of this bilinear form obviously has 
the following properties: every element below the main diagonal is zero; 
every element in the upper right hand fourth of the matrix is zero, the 
division being made by horizontal and vertical lines through the center; 
the determinant of the matrix has the non-zero value { l n (a) l n (b)\ n . 
Now this bilinear expression may be written in the normal form 

(6) [_R(u, v) + \B 1 (u, v)1zl = E U t (u) Vide), 

where the Ui and the V, are linear and homogeneous functions with con- 
stant coefficients-, the first in the variables u and the second in the variables 
v. From the totality of such representations we propose to choose those 
having certain properties. 

In case L\(u) = g u, so that our bilinear form is independent of X, 
the linear forms Z7< may be chosen in any way so that the 2n forms are 
linearly independent. The forms Vi are then uniquely determined and are 
linearly independent. This is the choice to be taken in the case of the 
Birkhoff problem. 

* If m = we take Ri(u, v) to be identically zero. 
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In the more general case when the bilinear form contains X we wish to 
realize the following conditions : the linear forms TJi, TJ 2 , • • • , U n , V„+i, 
V n+2 , • • • , Vi n shall be independent of X; the whole set TJ\, TJ 2 , • • • , TJ 2n 
shall be linearly independent for every value of X; the whole set V\, V 2 , 
• • • , V in shall be linearly independent for every value of X. For the last 
two conditions it is necessary and sufficient that the determinant of the 
set of linear forms TJi shall be independent of X (except possibly for a factor 
dependent on X but vanishing for no value of X) and different from zero, as 
one sees readily by aid of the fact that the determinant of the TJi multiplied 
by the determinant of the Vi has a value equal to the (non-zero) value of the 
determinant of the bilinear form. That these three conditions may in 
fact be realized may be seen from the following particular choice of linear 
forms: 

Ui = w<"-«(a), TJ 2 = u<"-»(o), • • •, U n -, = u'(a), U n = «<"-«(&); 
V n+ i = «<»-»(o), V n+2 = v'(b), V n+S = »"(&), • • ., V 2n = «<«-»(&); 

for i si n, Vi = total factor of TJi in the fundamental bilinear form; for 
i > n, TJi = total factor of Vi in the fundamental bilinear form. Indeed 
it is now evident that the needful properties may be realized in a great 
variety of ways. 

When the forms TJi and Vi in (6) are determined in accordance with 
the conditions named, then equation (1) and conditions (2) set a given 
boundary value problem the adjoint of which is set by equation (3) and 
conditions (4). It is clear, moreover, that problem (]), (2) is likewise the 
adjoint of problem (3), (4), so that the relation between the two problems is 
a reciprocal one. 

Remark 1. — We have seen that the conditions named can always be 
realized whatever the value of m, subject merely to the condition that it 
shall be less than n. For smaller values of m we also have certain addi- 
tional freedom. We may more generally take k linear forms TJi and 2n — k 
linear forms Vi independent of X where k is not necessarily equal to n but 
may be less or greater by an 'amount depending on m. When m is zero k 
may be any positive integer less than 2n* We may then make correspond- 
ing modifications in the two adjoint boundary problems. For the purposes 
of this paper, however, it is desirable to realize the maximum symmetry in 
the relations of the two problems; and consequently we confine our atten- 
tion to the case indicated in the definition as explicitly stated above. 

Remark 2. — It is proper to inquire whether the restriction that m shall 
be less than n is essential. The case m> n reduces to the case m < n 
by replacing X by 1/X (except for the particular value X = 0), so that we 

* See a treatment of this case by B6cher, Transactions of the American Mathematical 
Society, 14 (1913), 403-420. 



236 Caemichael: Boundary Value and Expansion Problems. 

need to consider only the case m = n. In this case the determinant of the 
bilinear form in the first member of (6) is not in general independent of X; 
and hence the properties of the forms Ui and Vi could not be maintained in 
general, properties which are needed in the further development. These 
properties might, however, be retained in large measure when m = n 
provided that g n and its derivatives have suitable values at a and at b; 
but we shall not attempt to carry this case along with the other, owing to 
the special hypotheses which would be necessary and the complications 
which would raise in certain parts of the work. 

Remark 3. — The question of generalization to an equation of the form 

L(u) + XZi(u) + X 2 £ 2 0) + • • • + \ k L k (u) = 

also arises. Under suitable hypotheses it is not difficult to maintain the 
necessary characteristics of the boundary conditions. Some of the theorems 
are still maintained without change;, but as early as theorem III certain 
difficulties would arise in general with this more extended case. 

We may now readily prove the following theorem : 

Theorem I. If for X = X a solution* u(x) of (1), (2) exists, then a solu- 
tion v(x) of (S), (4) also exists for X = X; if u(x) is unique (except for a con- 
stant factor) , v(x) is also unique (except for a constant factor) . 

Since u is a solution of (1), (2) for X = X, we have 

Ux(u) = U 2 (u) = • • • = U«(u) = 0, 

and for some k 

U n+k (u) =f= 0, 

where in U„+k we have the value X of X, since otherwise we should have 
u(a) = «'(o) = • • • = «<»-«(o) = u(b) = u'(b) = • • • = «<«-»(&) = 
and our function u(x) would be identically equal to zero, contrary to hy- 
pothesis. 

In the n-fold linear spread of solutions v(x) of (3) for X = X, there is at 
least one, say v(x), which satisfies the n — 1 linear homogeneous conditions 

V n+ i(v) = 0, 

where i runs over all the numbers of the set 1, 2, • • • , n except k. Now 
we have L(u) + XZi(w) = and M(?) + Xifi(») = 0. Hence if we mul- 
tiply the second identity in (5) by X and add to the first we find on in- 
tegrating that R(u, v) -j- \Ri(u, v) is a constant. Therefore from (6) we 
see that V ' n+k (u)V „+*(£) = 0; whence it follows that V n +k(v) = 0, and 

*By "solution" we usually mean a function which satisfies the equation and condi- 
tions but is not identically equal to zero; when we intend otherwise we shall indicate that 
fact explicitly. 
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hence that all the conditions (4) are satisfied by v, so that it is in fact a 
solution of (3), (4). 

It remains to be shown that v is unique whenever u is unique. For 
this purpose let us suppose, if possible, that v is not unique. Then two 
linearly independent solutions of (3), (4) exist, say v and v*. Then k and j 
exist such that 



V k (v*) Vj(v*) 



+ 0; 



for, otherwise, constants c and d would exist,* not both zero, such that 
Vi(cv-\- dv*) = for I = 1, 2, • • •, n, whence it would follow that the 
solution v = cv + dv* would satisfy all the conditions F,-(t>) = 0, i = 1, 2, 
• • • , 2n, and this is impossible since the V's are linearly independent for 
every value of X. Now choose u*, linearly independent of u, so as to satisfy 
the n — 2 conditions 

Ui(u*) = 

for i running over the numbers of the set 1, 2, • • • ,n, except k and j. Then 
from (5) and (6) we see that we have 

U k (u*)V k (v) + Uj(u*)Vj(v) = 0, 
U k (u*)V k (v*) + Uj(u*)Vj(v*) = 0. 

From this we see that U k (u*) = = Uj(u*). Hence w* satisfies (1), (2). 
Therefore, if the solution of (3), (4) is not unique for X = X neither is the 
solution of (1), (2) unique. Hence we conclude to the truth of the second 
part of the theorem. 

A value of X for which the system (1), (2) [[and hence the system (3), 
(4)3 has a solution will be called a characteristic value. The characteristic 
value is said to be simple if the solution corresponding to it is unique (except 
for a constant factor). 

Theorem II. If y\, ?/ 2 , • • •, y n are n linearly independent solutions of 
(1) for X = X, the condition that X is a characteristic value is that the deter- 



* In more detail one may proceed thus: distinct numbers hi and ji exist such that 

V h (7) 4= 0, V h (v*) * 0. 
If then 

V h (v)V h (v*) = V h (v*)V h (v) 
we have _ 

V h (v*) 4= 0, V h (v) 4= 0. 

But if our determinant is always zero it is sufficient to take 

c = V k .(v*), d = - V h (v), 
in order to realize the result stated. 



238 Carmichael: Boundary Value and Expansion Problems, 

minant 



A = 



tflfol) Ufa) ••• Utfjfn) 

U 2 (jfi) U 2 (y 2 ) ■■■ U 2 (y n ) 
U n (yi) U n (y 2 ) ••• Un(y n ) 



shall vanish; a characteristic value X of X is simple when and only when some 
first minor of A does not vanish. 

If we take the general solution of (1) in the form 

u = ayi + c 2 y 2 + • • • + c n y n , 

we see that the vanishing of A is a sufficient condition that conditions (2) 
may be satisfied through an appropriate choice of c\, c 2 , ■ •-, c„. When 
some first minor does not vanish this choice is unique (except for a constant 
factor). 

Theorem III. If u t {x) and v 3 {x) are solutions of (1), (#) and (3), (4), 
respectively, the first for X = X, and the second for X = Xy, and if X* 4= Xy, 
then we have 

Jl b UiMi(vj)dx = 0, f£vjL\(Ui)dx = 0. 

From the equations 

L(ui) + \iUiUi) = 0, M(vj) + XiJfi(*i) = 
we have 

{v£L(Ui) + \iI*(Ui)l - ulM(Vj) + \iMi(Vj)~]} + (Xi - Xy)«iMi(«y) = 0. 

If we integrate from a to b we see through aid of (2), (4), (6) that the expres- 
sion inclosed in braces yields the value zero. Hence we have 

(Xi - \j)f a b UiMi(Vj)dx = 0. 

Dividing through by the factor X, — Xy we have the first relation in the 
theorem. The second may be proved in a similar way. Moreover, either 
of these relations follows readily from the other by means of the second 
relation in (5) and relation (6), looked upon as an identity in X, together 
with the boundary conditions. 

If we have two sets of functions Ui(x), u 2 (x), u 3 (x), • • • and Vi(x), v 2 (x), 
v$(%), •••> and the associated adjoint linear forms Li(u) and Mi(v) such 
that relations (7) are satisfied for every i and j which are different we shall 
say that the two sets of functions are adjoint each to the other. If the two 
sets are the same and the form L\{u) is self -ad joint we shall say that the 
single set is self-adjoint. The conceptions thus introduced are obviously 
generalizations of the conceptions of biorthogonal sets and orthogonal 
sets of functions. 

If now we have further 

f a b u i M 1 (v i )dx #0, t = 1, 2, 3, • ■ •, 
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it is clear that we may formally determine the coefficients c,- in the formal 
expansion of a given function f(x) in the form 

f(x) = CiUi(x) + c 2 u 2 (x) + •■-. 

In fact, if one multiplies by M^vi) and integrates from a to b one has for- 
mally 

i= 1,2,3, •••. 



c _ fa b f(x)M 1 (.V i )dx 

f a b u i M l {io l )dx 



When X is not a characteristic number for (1), (2) there exists* a unique 
0(x, s; X) such that the solution <p of 

L(<p) + XZi(p) = r, CTi(p) = U,(<p) = • • • = *7»(«») = 0, 

is given by 

<P = fa b G{x, s; \)r(s)ds; 

and likewise there exists a unique H(x, s; X) such that the solution ip of 

1(f) + XilfiCf) = r, F^-iW = *Wtf) ='" = Mf) = 0, 

is given bv 

f = /„ 6 #(a:, s; X)r(*)d». 
Moreover, we have 

G(x, s; X) = #(s, x; X). 

Explicit formulae for G and # are readily obtained (see Bocher, I.e.). 
The formula for G may be written in terms of a fundamental system y u y 2 . 
■ • •, y n of solutions of (1) as follows: 

yi(x) y 2 (x) • • • y„(x) G(x,j; X) 
U 1 (y 1 ) U,{y 2 ) ••• U,{y n ) U,{G) 



(8) G(«, s; X) = 



N(x, s; X) 



Un(yi) U n (y 2 ) ■■■ U n (y n ) U n (G) 



A(X) A(X) 

where A(X) has the value given in theorem II and G(x, s; X) has the value 



(9) ■Q(x,s;X)= igj- 1 



2/i(a;) 2/ 2 (a;) 
2/i ( "- 2) («) 2/2 ( "- 2) (5) 



2/l(«) 2/2 (*) 



2/»(a0 
2/« (n - 2) («) 

Vn(s) 



(*) 



2/i ( " _1) («) 2/2 (n_1) (s) 
2/i (n_2) (5) 2/2 (n - 2) (5) 



y. M («) 



2/n 



(«-2) 



2/lW 2/ 2 («) ••• 2/n(«) 

the sign + being taken when x > s and the sign — when x < s. 
* B6cher, "Les Methodes de Sturm," Chapter V. 
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Now the functions y\, y 2 , ■ ■ • , y n may be taken analytic in X. Hence 
G is analytic in X except when A(X) = 0, that is, except for the characteristic 
values X, of X for problem (1), (2). If X; is a simple characteristic value 
we may write 

(10) G(x, s; X) = ?^±+ a {x, s; X) 

A — A; 

where a{x, s; X) is analytic in X for X = X*. From (8) we have 

„ , , _ N(x, s; Xj) 
H t {x, s) - ^^ 

where A'(Xi) denotes the value at X = X» of the derivative with respect to X 
of A(X). In N(x, s; X») the coefficient of G is zero. Hence N(x, s; X,) 
and its first n derivatives with respect to x are continuous. Moreover, it is 
easy to see that it satisfies (1), (2); therefore Ri(x, s) is a solution of (1), 
(2). From the relation between G and H we infer that Ri(x, s) is a solution 
of (3), (4) in s. From the definition of simple characteristic value it follows 
then that 

Ri(x, s) = CiUi(x)Vi(s) 

where c» is independent of x and s and is yet to be determined. 
Now from (10) we have 

lim {(X — \i)G(x, s; X) — c,-w,-(a;)fl»(*) } = 0; 

whence it follows that 

(11) lim [(X - \i)f a h G{x, s; \)I*(in(s))d8 - c l Mx)f a %(s)L 1 (u i (s))ds'] = 0. 

A=A» 

But since 

L(u f ) + \Li(ui) = (X — Xi)Zn(Mi), Ui(ui) = U 2 (ui) = • • • = U n (ui) =0, 
we have from the fundamental property of G that 

(X — 'Ki)J > a b G(x, s; \)Li(Ui(s))ds = Ui(x). 

Substituting this value into equation (11) we have 

CiJl h i}i{s)Li{Ui{s))ds = 1. 

Similarly we should obtain the relation 

dfJ'u^MiiViis^ds = 1. 

Hence we have for the residue Ri(x, s) the value 

Ui{x)Vi(s) _ Ui(x)Vi(s) 



(12) Ri(x, s) = 



If r is a contour in the X-plane which incloses (just) the characteristic 
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numbers Xi, X 2 , • • • , X n , and if the latter are simple, we have readily 
■■ fvfa b G(x, s; X) , , f(s)dsd\ 



K ' = f> fa b f(s)M 1 (n(s))ds 

fet f a "u i (s)M 1 (v i (s))ds A) ' 

Thus we have represented in the form of a contour integral the sum of 
a finite number of terms of the formal expansion 

fci f a h Ui{s)Mi{xii{s))ds 

of a function /(a;) for the case when the characteristic numbers Xi, X2, X3, 
• • • are all simple and G(x, s; X) has a pole of the first order. At a value 
X; of X for which these conditions are not satisfied, the corresponding term 
in the formal expansion in (14) should be replaced by 

where Ri(x, s) denotes the residue of G(x, s; X) at X = X* and Vi(s) is the 
solution N(x, s; X,) of (3), (4) in s. Thus we have a definite formal expan- 
sion for every given function f(x) for which the integrals involved exist. 

In accordance with the general plan of this paper we rest this investiga- 
tion at this point. We have carried it far enough to set clearly the funda- 
mental expansion problem associated with it, namely, the problem of deter- 
mining those functions f(x) for which the given formal expansions are valid. 
This problem Birkhoff (I.e.) has treated with great success for the special 
case considered in his paper. It appears likely that the methods of 
Birkhoff's paper are sufficient to carry through the more general investiga- 
tion. The results obtained above are sufficient to initiate the further 
treatment. 

In the following sections we derive formal expansions in connection with 
a great variety of boundary value and expansion problems. We do not 
go into as much detail as in this section; but the reader who has this section 
in mind will have no difficulty in seeing how the results actually obtained 
in each case set the fundamental expansion problem for that case. 

2. The Problem Arising from a System of Differential Equations of the 
First Order with One Parameter and from Related Integro-differential 
Equations. — Let us consider the adjoint systems of differential equations 

dv ■ n 

(1) dx = ^ ( aij + Xa ^'> i = 1, 2, • • •, ra, 

dz- v^ 

(2) -j 1 = ]T (— an — \aji)z 3 ; i= 1, 2, • • • , n. 
ax fel 
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If for fixed i we multiply these respectively by z» and yt, add the resulting 
equations member by member, sum as to i from 1 to n, and then integrate 
from a to b, we have 

(3) \j/izi + y 2 z 2 + •■• + 2/»z«]S* = 0. 

Let us now suppose that homogeneous linear boundary conditions 
(involving the end-points a and b) are set up on the y's and adjoint condi- 
tions on the z's so that (3), for the solutions to be considered, is satisfied in 
virtue of the boundary conditions alone. Moreover, let us suppose that this 
is done in such a way that the characteristic values of X for (1) with its 
conditions are the same as those for (2) with its conditions, the situation 
being analogous to that treated in § 1, and that the number of characteristic 
values is infinite. Let us denote the distinct characteristic values and 
corresponding solutions by the symbols X w , y/ k) , z, (i) , k = 1, 2, 3, • • •• 
In (1) let us now replace X and y { by X w and y^; in (2) let us replace X 
and Zi by X (0 and zP, where X (0 #= X w ; then we have two new systems. 
For fixed i let us multiply the former by z» (0 and the latter by yf- k) , add, 
sum as to i from 1 to n, integrate from a to b, and simplify the first member 
by means of the equation to which (3) reduces when y t and z t are replaced 
by yi m and z» (0 . Thus we have 

= fa b il ±, [><,■ + X<*>a<,)yy<»z«<«> - (o,< + ^a^y^z^dx, k 4= I; 

whence, on simplifying and dividing by X (i) — X (0 , we have 

(4) Sa b t, t, a jiy ^W l) dx =0, k +Z . 

If we suppose that the integral in the first member of (4) is different in 
value from zero when k = I, we have a convenient means for the formal 
determination of the coefficients c k (independent of x and i) in the simul- 
taneous formal expansions of n given f unctions /i (x), / 2 (x), • • -,/«(«): 

(5) /<(*) = E c^» a) , i = 1, 2, • • •, n. 

For determining these constants c k we proceed as follows: multiply both 
members of (5) by ay#/°, sum as to i and as to j from 1 to n, and integrate 
from a to b; in this way we obtain a relation from which we have readily 
the value of ci, namely: 

n n 

f a h H J^ajifiZj^dx 

(6) c ; = P-^ 1=1,2,3, .... 

A b Y, H^yPz^dx 
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Thus we have in (6) the formulae for the coefficients c*, the same for each 
function fi(x), needed in the joint expansions in (5) of the n functions fi(x). 
They are formula? of great interest. 

We have outlined very briefly the foregoing expansion problem, which 
is probably the same as that of an unpublished paper of Birkhoff (see 
abstract in Bulletin of the American Mathematical Society, (2) 25 (1919), 
p. 442), because we wish to utilize the suggestion afforded by it for setting 
up a related problem for certain integro-differential systems obtained from 
the foregoing differential systems by the limiting process which Volterra 
has so frequently employed and to such great advantage. It is clear that 
the notion of joint expansion utilized in connection with equation (5) is an 
essential element in the formulation of the problem. This notion I obtained 
not from Birkhoff's abstract (which is too brief for details of any sort) but 
from the dissertation of Dr. C. C. Camp which the author was kind enough 
to allow me to read in manuscript.* In this dissertation there is a detailed 
and full treatment of a very special case of the expansion problem formulated 
above; and the range of ideas developed has been of great use to me in 
the formulation of the boundary value and expansion problems of this paper. 

If we apply to systems (1) and (2) the Volterra limiting process we may 
obtain the integro-differential equations 

(7) ^^ = fj {u(x, s, t) + Mx, s, t) }y(x, t)dt, 

ox 

(8) ^-^ = J? { - u(x, t, s) - \v(x, t, s) }z(x, t)dt, 

ox 

where the range of variation of x is from a to b while that of s and t is from 
a to (8. Let us multiply both members of (7) by z(x, s) and both members 
of (8) by y (x, s), add the two resulting equations member by member, and 
integrate with respect to s from a to /?; then, on carrying out certain 
simplifications through an interchange of the order of integration with 
respect to s and t and through certain obvious reductions, we have formally 
the relation 

fj ^ {y(x, s)z(x, s) }ds = 0. 

Integrating with respect to x from a to b, we have 

(9) f/{y(b, s)z(b, s) - y(a, s)z(a, s)}ds = 0, 

a relation which corresponds to (3) for the differential systems. 

* Since this was written a paper by A. Schur has appeared dealing with this problem; 
see Mathematische Annalen 82 (1921) : 213-236. It contains (p. 216) a reference to a 
special case of the problem treated by Hilbert (Gottingen Nachrichlen, 1906, pp. 474-480) . 
This I had not seen before. 

Dr. Camp's article will appear in the next volume of this Journal. [Editor.] 
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The first member of the last equation is a sort of bilinear form in a non- 
enumerably infinite set of variables y{a, s), y(b, s), z(a, s), z(b, s), a set of four 
variables being formed for each s of the interval (a/3). One desires to have 
a reduction of this to a normal form analogous to that used in connection 
with equation (6) of § 1. At present we do not attempt to solve this 
problem in general. What we desire is to have a set of boundary conditions 
on y and z separately, linear and homogeneous in character, so that these 
boundary conditions alone shall imply equation (9) and so that the y- 
problem and the z-problem shall have the same characteristic values X 
infinite in number. Assuming that such equations have been set up, we 
seek to determine the formal character of the boundary value problem in 
its simplest aspects. 

Let us represent characteristic numbers and corresponding solutions 
by Xjt, yk(x, $), Zk(x, s), k = 1, 2, 3, • • • . Let us replace A and y in (7) by 
Xfc and yjc', let us replace X and z in (8) by X; and zi, k =|= I; let us multiply 
the first resulting equation member by member by zi(x, s) and the second 
by yk(x, s), add the two resulting equations member, integrate with respect 
to s from a to (3 and with respect to x from a to b, and simplify the resulting 
relation by aid of what equation (9) becomes when y and z are replaced by 
yic and zi respectively. From both members of the equation which so 
results remove the non-zero factor X& — X;; thus we have 

(10) JVS/SM*, t, s)y h (x, s)zi(x, t)dsdtdx =0, k * I. 

If we suppose that the first member of (10) is different from zero for 
k = I and if we seek constants cj- so that a given function f(x, s) shall have 
the formal expansion 

(11) fix, s) = X) c k y k ix, s) 

k=\ 

we are led (through the relation obtained on multiplying both members of 
(11) by vix, t, s)zi(x, t) and on integrating appropriately) to the formulae 
for Ci, namely: 

/ 12 ) Cl = flf^fl^! t> s )ffa s)zi{x, t)dsdtdx t 1= 1 2 3 ■■■ 
Jl b fJ/'J''>ix,t,s)yi(x,s)ziix,t)dsdtdx 

As a special case of the foregoing problem let us take that in which u 
and v are independent of x and let us write 

u{x, s, t) = g(s, t), vix, s, t) = his, t). 

Solutions exist of the forms 

yix, s) = pix)<r(s), zix, s) = R(x)Sis), 
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where each function in the second members is a function of the single 
variable indicated. As special boundary conditions implying equation 
(9) we may take the following: 

y(b, s) - y(a, s) = 0, z(b, s) - z(a, s) = 0. 

These conditions now reduce to the following : 

p(b) - p(a) = 0, R(b) - R(a) = 0. 

If we substitute the foregoing special values of y and z into the special 
forms of (7) and (8) indicated we may separate variables in the way which 
is classic in the study of the partial differential equations of physics and so 
come through to the following relations : 

, . p'(x) = pp(x), p(b) - p(a) = 0, 

^ ' P<r(s) = fj {g(s, t) + \h(s, t)}a(t)dt; 

R'(x) = - pR(x), R(b) - R(a) =■- 0, 
{ > - pS(s) = fJ \9(s, t) + XA(*, t)Wt)dl, 

where p, is a constant. In (13) and (14) we use the related constants p 
and — p so as to maintain the requisite property of adjointness. The 
differential system in each of these pairs of conditions is easily solved, 
appropriate characteristic values of p being thus determined. Then the 
integral equation of the pair is to be solved, appropriate characteristic 
values of X having been determined for each value of p. 

Let us further specialize. Suppose that we take b = w and a = — ir. 
Then for p we have the values raV— 1, where n is a positive or negative 
integer, and for p(x) and R{x) the corresponding values e nxV ~ l and e _M; ''-i 
respectively. If we take g(s, — and employ a non-zaro value of p 
we have in the second line of (13) and in the second line of (14) an integral 
equation with parameter X/V in the one case and — X/ju in the other case. 
It is well-known that the quantities a, jS, h(s, t) in this equation may be 
so chosen that for each value of p corresponding values of X exist giving 
rise to the functions cos ms for solutions where m runs over the set 1, 2, 3, 
• • •. For the special case thus defined the expansion (11) has the form 

00 

f(x, s) = 2 ( c »» cos mse nix + y mn cos mse - "**). 

m, n—\ 

If f(x, s) is a real-valued function when x and s are real it is clear that 
c m „ and y mn must be conjugate imaginaries. Then the expansion reduces 
to the form 

oo 

f(x, s) = 2Z ( a mn cos nx + b mn sin nx) cos ms. 
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where a mn and b mn are real numbers. Furthermore, if f(x, s) is an even 
function of x as well as of s we must have b mn = so that the expansion 
becomes the double cosine series 

CO 

f(x, s) = 2^ o n cos ms cos nx. 

Thus we see that the expansions afforded by (11) and (12) contain as special 
cases some of the classical expansions of analysis. 

The problems formulated in connection with (1), (2) and (7), (8) are 
capable of various generalizations. It will be sufficient to indicate briefly 
one of them. Let us consider the system of equations 



(15) 



^^=t,(a ij +\a ij )y j (x > s) 
ox j=l 



+ f/JXpa(x, s, t) + \<7ij(x, s, t)~]yj(x, t)dt, 



3=1 



dzi(x, s) 



(16) °^lll = £ (- a H - \a H ) zj(x, s) 

ax j=i 

j=i 

f or i = 1, 2, • • • ,n. Here in general the a# and a„ are functions of x and 
s. In the special case when p# = = o-*y and a tj and a i} - are functions of x 
alone, we may take yt and 2, to be functions of x alone, and we have our 
systems (1) and (2). It is clear that equations (7) and (8) are included in 
another way as a special case. 

Let us multiply (15) by z { (x, s) and (16) by yi{x, s), integrate the resulting 
equations as to s from a to j8 and in the relation coming from (16) interchange 
the order of double integration and in the part of the formula containing 
this double integration interchange s and t, add the two equations thus 
gotten from (15) and (16), and in the resulting equation sum as to i from 
1 to n. Finally in this resulting equation integrate as to x from a to b. 
Then we have the relation 

(17) SSEtiViQ, »)*&, «) - V*(a> s >i( a > *) }> = 0- 

This corresponds to our previous relation (9). 

Let us suppose that linear homogeneous boundary conditions on the 
?/; and the z 4 - separately are so chosen that relation (17) is implied by these 
boundary conditions, that the ^/-problem and the 2-problem have the same 
characteristic values X infinite in number, and let us denote the correspond- 
ing values and solutions by X&, 2A W , z» (;fc) , k = 1, 2, 3, • • •. Then by a 
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method of procedure which is now obvious we come through to the funda- 
mental reciprocal relation 

fa b {S/rit EM*, *)yi w {x, *)*i m (x, «» 

(18) i=lj =\ n 

+ S/S/CZ 2>*(*, t, »)vi w (x, 'WHx, t)2dsdt\dx = o, k*L 

This is the generalized biorthogonality condition of which we have already 
found special cases. 

Let us denote by D k i the first member of (18) and let us suppose that 
Dkk 4= for each k. Then if we have a set of expansions of the form 

(19) fi(x, s) = I>r2/i (r) fo s), i= 1, 2, •■■,», 

r=] 

giving the simultaneous expansions of the n functions fi(x, s), i = 1, 2, 
• • •, n, in terms of the y^'ix, s), the coefficients c r may be found by mul- 
tiplying both sides of the equation 

n n 

£ a H (x, sWKx, s) + f/CZvjiix, t, sWHx, t)yt, 

summing as to i from 1 to n, integrating as to s from a to j8 and as to a; from 
a to 6, and simplifying by use of (18) ; thus we have 

V vv c v = Sa b {S/Ct £,<**(*, »)M*. »WK*, *) 

(20) . . fcl *= x 

+ SJT, 2>w(*» t, *WKx, t)fi(x, s)dt]ds}dx, v = 1, 2, 3, • • •. 

It is clear that one may apply to (15) and (16) the Volterra limiting 
process and so replace these equations by others in which the discrete 
variable i has given way to an additional continuous variable, the resulting 
integro-differential equation containing a repeated integration, and that the 
formal expansion problem subsists in the new form. One might then employ 
a system of such equations and apply again to them the Volterra limiting 
process; and so on indefinitely. 

3. Differential and Integro-differential Problems Involving v Param- 
eters. — Let us consider the v systems of equations 

dii ■ "* 

,, . — j— = Z-i (dOkij + Xlffllitf + \2a2kij + - - - ~f" ^v a vkii)ykh 

(1) dx i=l 

i = 1, 2, • ••, raj,, 

a separate system being formed for each value h of the set 1, 2, 3, • • •, v; 
and with these let us associate the v systems of adjoint equations 

dz • nk 

(2) — r— • = 2Li ( — %)< — ^i a ikji — • • • — \ v a v kn)zkj, i = 1, 2, • • •, ?ifc. 

eta* .?=i 
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For fixed k and i multiply equation (1) member by member by z ki and equa- 
tion (2) member by member by vm, add the two resulting equations member 
by member, sum as to i from 1 to n*,, and integrate with respect to x from 
a to b; thus we have 

(3) ZyklZkl + Vk&kl + ■ • ■ + Vkn^kn^zl = 0, k = 1, 2, • • • , V. 

Let us now suppose that for each value of k we set up on the quantities 
Vki(b), ykiifl), i = 1> 2, ■•-, n k , and separately on the quantities Zki(b), 
Zki(a), i = 1, 2, • • • ,n k , linear homogeneous boundary conditions (analogous 
to those developed in detail in § 1) of such sort that relations (3) are implied 
by the boundary conditions alone and that the ^-problem and the z-problcn 
have the same sets Xi w , X 2 (s) , • • •, X„ w of characteristic values, that the 
number of sets is enumerably infinite and that these distinct sets are 
denoted by the named symbols for s = 1, 2, 3, • • • . Let the corresponding 
solutions be denoted by yki (s) , Zki W - 

In equation (1) let us replace the quantities y and X by corresponding 
quantities y M and X (s) ; in equation (2) let us replace the quantities z and X 
by corresponding quantities 2 (r) and X (r) where r #= s. For a fixed k multiply 
the first resulting equation member by member by zm w and the second 
by yu (s) , add the two resulting equations member by member, sum as to i 
from 1 to n k , integrate with respect to x from a to b, and simplify by what 
(3) becomes on replacing the quantities y and z by the corresponding quan- 
tities y M and z (r) . In this way we are led to the following relations: 

(4) (XxW - Xi«)y; 6 Z EoiwWW'Ae + • • • 

i=l ,7=1 

m* nit 

+ (X,« - K M )Sa b Z Za vm yki M z kj M dx = 0, 

4=1 .7=1 

a separate relation being formed for each value k of the set 1, 2, • • •, v. 
Now the differences Xi (s) — Xi (r) , • • •, X„ (s) — X„ (r) are not simultaneously 
zero; hence the determinant of their coefficients in the system (4) must 
have the value zero, that is, we must have 



Jl h Tl zL a ikjiyki (s) Zk} ( - r) dx 



= 0, r * s, 



where the term of the determinant explicitly written is that in the Mi 
column and fcthrow. 

The last relation may be reduced to a more convenient form. In the 
kth row of the determinant let us replace the variable x of integration by 
tk and the variables i and j of summation by i k and jk respectively. Then 
the relation reduces to 



rv rb, . . ri 

J a J a Ja 






dtidti • • • dt v = 0, r #= s, 
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the determinant employed being that whose element in Zth column and kih 
row is that explicitly written. Then if we employ the symbol D to represent 
the determinant 

(5) Di 1 j 1 ... ivj r (ti, t 2 , • • • , t v ) = \aikj k i k (h) I 

whose element in Zth column and &th row is that explicitly written, we have 
for the final form of our generalized biorthogonality conditions the following: 



(6) 



2-i 2-i ' " ' 2—i 2-i "i\h ••• w»w v 



■,ty) 



h=L J 



r =(= s. 



In what follows we shall suppose that we have before us the case in which 
the first member of (6) is different from zero when r = s. 
Now let us consider the set of functions 

Ui, - iy(h, h, •••, t v ), i k = 1, 2, • • •, n k for k= 1, 2, ■■-, v, 

with regard to the question of their simultaneous expansion in the form 

(7) Ai a - iy(tl, U, ■ ■ -, ty) = £C4^(*1)^««) • • • M*>)> 

8=1 

the coefficients c s being the same for every function of the set. In case such 
an expansion exists we have a ready means of determining the coefficients 
c s . The form of the first member of (6) suggests the operations to be 
performed. Let us denote by B rs the first member of (6) and let us sup- 
pose that B rr #= for r = 1, 2, 3, ••• . Then for the determination of the 
c r we have the relations 



r Ml m 
h rr C r = J a J a • • ' J a \ / . 2—1 . 



(8) • • •2J2JDi 1 j 1 ••• iyjy(h, h, • • • ,K)fhH ••• iy(h, h, • •• , ty) 

iy=lj!=l 

■$l{h)z%M ■ ■ ■ sSU0}*i<& • • • dty, r = 1, 2, 3, ■ ■ ■ . 

The expansion problem which is set by these formal results we may 
look upon as a generalization of that briefly treated at the beginning of § 2. 
From the fuller treatment of a related problem given in § 1 one sees readily 
how to proceed to a more precise formulation of the problem here briefly 
characterized. Moreover, it is clear that one may formulate the problem 
at once for a system of v equations of general orders and so include in one 
form (rather complicated, to be sure) the three differential equations 
problems already set. 

We saw in § 12 how the integro-differential equations problem (7), (8) 
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grew out of the differential equations problem (1), (2). In a similar way 
the problem (1), (2) of this section gives rise to a new problem associated 
with the equations 

dy k (x, s) = „ ^ + Xi0i s>t)+ ... 

(9) ox 

+ \ v a vk (x, s, t)}y k (x, t)dt, 

/im — iH — = f* i ~ a ° k ( x > f > s ) ~" A i a xit(a;, *, s) 

(10) ox 

— • • • — \a yk {x, t, s)}z k {x, f)dt. 

With the guide afforded by the earlier part of this section and the special 
case of our present problem treated in § 2 one may proceed readily to set 
up the generalized biorthogonality conditions and by aid of them to deter- 
mine the coefficients in the formal expansion of a function of 2v variables. 
Corresponding to equation (6) we should thus have 

fa h fa h - ■■Sa b {f/S/ ' ' -f/D(x U X 2 , • ■ • , X„ » U h, »%, U, • • •, *„, O 

• II ^(^ sa)za°(«a> h)dsidh • ■ • ds v dt v \dxidx 2 ■ • ■ dx v = 0, r 4= s, 

where the symbol D denotes the determinant of order v whose element in 
7th column and &th row is ai k (x k , s k , t k ). If we denote the first member 
of this relation by C rs and suppose that C rr 4= for every r, then for the 
coefficients c r in the expansion 

00 V 

(12) f(xi, xi, ••-, x v , si, s 2 , • • • , s v ) = z^c.JJ yh (s) (xh, $h) 
we have formally 

O rr C r = J a J a ' ' ' Ja \Ja J a ' ' ' J a. *s\%l> ' ' 't X„, S\, l\, • • •, S„, t v ) 

v 

(13) /(.-ci, • • •, x v , si, • • •, s v ) -Ylzh\x h , t h )-dsidti • • • ds„dt„\ 

X dxidxi • • • dx v , r = 1, 2, 3, 

These results might, if it were desired, be also extended by means of the 
suggestion afforded by the problem associated with equations (15) and (16) 
of § 2 and the generalizations of it which are mentioned there. 

4. Problems Involving Difference and Integro-difference Equations. — 
Let us consider the adjoint systems of difference equations 

n 

(1) Ui(x + 1) = 2(5,7 + <pn + \\pij)uj(x), i = 1, 2, • • •, n, 

j=\ 

n 

(2) Vi(x) = £(5y.- + <f>ji + \^n)vj(se +1), i = 1, 2, • • •, n, 
where 5„ has the value 1 or according as i is or is not equal to j and <py 
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and \pij are functions of x which are analytic at infinity and vanish there 
to an order at least as high as the second. Under such hypotheses these 
equations have for a fixed value of X fundamental systems* 

uij{x), u 2 j(x), • • •, u nj (x); vij(x), v 2 j(x), ■••, v nj (x) 

of solutions each function of which is analytic throughout the finite plane 
except for a set of singularities in a left half -plane; moreover, if x approaches 
infinity along any ray from the origin exclusive of the negative axis of reals 
or along any line proceeding to the right parallel to the axis of reals, we have 

lim Uij{x) = 5,y, lim «#(*) = 8#. 

Any solution of either equation formed from the foregoing solutions of that 
equation by taking a linear combination of such solutions with constant 
coefficients will have the property that each function in it approaches a 
limiting value as x approaches infinity in the manner specified. We confine 
attention to such solutions of equations (1) and (2). 

Now if we multiply equation (1) member by member by Vi(x + 1) and 

(2) by — Ui(x), add the resulting equations member by member, and sum 
as to i from 1 to n, we have 

(3) EA|«.W»iW! = 0. 

If the real part of a is sufficiently large we have u t and «,■ analytic at every 
point x whose real part is not less than the real part of a. Hence in (3) 
we may sum as to x from a to infinity, where x runs over the values a, a + 1, 
a + 2, • • • ; thus we have 

TO 

(4) £{«;(<*> )»;(<*>) — Ui(a)Vi(a)} = 0. 

This is clearly analogous to equation (3) of § 2 and differs from it (apart 
from notation) only in having the limit co instead of the limit b. 

Let us now suppose that adjoint homogeneous linear boundary condi- 
tions are set up on the w»(°o), w,(a) on the one hand and on the t>»(°o), »i(a) 
on the other hand so that relation (4) is implied by the boundary con- 
ditions and that this is done in such a way that the characteristic values 
of X for (1) and the conditions on the Ui are the same as those for (2) and 
the conditions on the Vi and that the number of characteristic values is 
infinite. Let us denote the distinct characteristic values and corresponding 
solutions by X w , Ui (k \ v/®, k = 1, 2, 3, 

In (1) let us now replace X and Ui by X w and w; w , in (2) let us replace 
X and Vi by X (0 and »» (0 , where I 4= k; then we have two new systems. Let 

* See American Journal of Mathematics, 35 (1913), 161-162. 
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us multiply the first resulting equation member by member by »/°(a; +1) 
and the second by — Ui W (x), add the resulting equations member by mem- 
ber, sum as to i from 1 to n, sum as to x over the set a, a + 1, a + 2, • • • , 
simplify by aid of what (4) becomes when u t and «< are replaced by u t w 
and »i (!) respectively, and divide by the non-zero quantity X (A:) — X a) . 
Thus we have 

(5) £ £ I>/.-(o + 0«* w (o + 0»/°(a +1 + = 0, k*l 

t=0 j=l j=l 

If we suppose that the quantity in the first member of (5) is (for all k) 
different from zero when I is replaced by k we shall have a ready means for 
the formal determination of the coefficients Cu in the (assumed) expansions 



(6) /<(*) = £ <*«,•<»(*), i=l,2, 



ji, 



the same coefficients Ck being employed for each of the n functions. For 
this purpose we multiply both sides of equation (6) by \f/ji{x)vj W (x + 1), 
sum as to i and j each from 1 to n, and sum as to x over the set a, a + 1, 
a + 2, • • • ; employing (5) we come through readily to the relations 



(7) 



£ £ £ **(« + 0/i(« + 0»/ B (o +1 + 
£ £ £ Ma + 0«* (B (« + 0»/ B (« + 1 + 



k = 1, 2, 3, 



Since the coefficients c^ are determined merely by the function-values 
at a discrete set of values of x it may seem at first sight that this is not 
likely to lead to an expansion theory of much interest. But it would be 
natural to confine attention to functions /»•(«) having suitably restricted 
properties, especially in the neighborhood of infinity, and for such functions 
it may be anticipated that the theory is valuable. Moreover, it furnishes 
also a natural means of interpolation of a function defined initially only 
over the set of points a, a + 1, a + 2, • • •, the function values at these 
points being the only ones needed in the determination of the coefficients 
Ck- Again, it affords a natural means of formal expansion of functions 
defined only for positive integral values of the argument, such as is the 
case with the usual number-theoretic functions. 

Let us rapidly illustrate the foregoing theory by considering the problem 
associated with the equations 

(8) «(3+l)-(l-^W) = 0, 

(9) (l-^W+1)- «(*) = (), 
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where /x 2 is the parameter. These equations have the solutions 

(10) u(x) = -i — ' \ '- , v(x) = — — \'f , 

lT(x)J T(x- [x)T(x+ m) 

where T(x) denotes the usual gamma function; and the general solution 
for each case is gotten by multiplying the particular solution by an arbitrary 
periodic function of period unity. But if u(x)[v(x)~] is to approach a 
definite limiting value when x approaches infinity in the sense defined 
above this periodic multiplier must reduce to a constant; and such constant 
multiplier is clearly irrelevant to the problem in hand. Hence, for our 
purposes, we take the solutions given in (10). 

Now condition (4) reduces in the present case to u(k>)v(<x>) — u{a)v{a) 
= 0. As suitable boundary conditions implying this relation we may take 

(11) tu(<x>) - u(a) = 0, jv(«>) - v{a) = 0, 

where t is a given non-zero constant. Now from (10) and the asymptotic 
properties of the gamma function we see that tt(oo) = 1 = jj(oo). Then 
the condition of consistency for the ^-problem and that for the ^-problem 
both reduce to the following : 

(12) T(a- v)T(a + n) = t{r(a)}*. 

The roots of this equation give the common characteristic values of n 2 
for the M-problem and the ^-problem. [For each solution ju of (12) there 
is also the solution — n; but only their common square ju 2 is a characteristic 
value for our problem.] 

For varying choice of a and t the expansion problem which results 
brings a variety of properties to notice. One of the most interesting cases 
is that in which a = \, — 1 ^ 1/t ^ -f- 1; in this case equation (12) 
readily becomes 

(13) cos w ix = - , 

if one reduces by aid of the relation r(*)r(l — s) = r/smws. Then the 
characteristic values ix 2 are readily determined in terms of a single one of 
them. Let ju be the smallest positive value of ix satisfying (13); then the 
characteristic values /x 2 are 

M 2 = (m + 2n) 2 , n = 0, ± 1, ± 2, ± 3, 

Corresponding to these let us denote the solutions by u n (x), v„ix) so that 
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we have* 

..... , . T(x - m - 2n)T(x + m + 2n) ,, 1 

(14) u n {x) = -i £ — =- > vv ~> »»(*) = — — , 

n = 0, 1, 2, 
The biorthogonality conditions now reduce to 

(15) «<*+*)'«.<*+ f)~ 0, + 

Now we may apply to systems (1) and (2) a Volterra limiting process 
analogous to that by which we arrived at equations (7) and (8) in § 2. 
By so doing we are led to the integro-difference equations f 

(16) Au(x, s) = S*{(p{x, s, t) + \\p(x, s, t))u(x, t)dt, 

(17) - Av(x, s) = JT{<p(x, t, s) + Mix, t, s) }v(x + 1, t)dt. 
Corresponding to equations (4) to (7) we now have the following: 

(18) f£ { u ( °° » S M °° » s ) ~ " u ( a > s)v(a, s)}ds = 0, 

(19) ifvMia+r, t,8)u k (a+r,»Ma+l + r,t)d8dt=0, k + I, 

(20) f(x, s) = J2 c k u k ix, s), 

k=\ 

E SSSMa + r, t, s)fia + r, s)v k (a + 1 + r, t)dsdt 

(21) c k = *-? , 

Z) /«V«V(°'+ r > t> s)u k ia +r, s)vkia + 1 + r, t)dsdt 

k= 1, 2, 3, 
The foregoing work of this section is analogous to a part of that of § 2. 
It is clear that the analogues of the other parts of § 2 exist here and can be 
readily developed. In fact, we can develop the theory of difference and 
integro-difference equations in forms analogous to all parts of the problems 
treated or suggested in §§ 1 to 3 for differential and integro-differential 
equations. 

5. Mixed Problems Involving One or More Parameters. — In view of the 
basic algebraic theory and the transcendental problems treated in §§ 1-4, 
it is clear that there must exist in the theory of integral equations expansion 
problems involving n parameters not only in the classic case when n = 1 
but also in the general case when n is any positive integer. It was my 
intention to exhibit the basic formulae on which such an expansion theory 
may be based; but, while this paper was in course of preparation, an abstract 

* It should be observed that the negative values of n merely give us the same solutions 
over again. It is also worth noting that u n (x)/ui(x) is a rational function of x so that from 
the expansion of f(x) in terms of u„(x) we have that of f(x)/v,i(x) in terms of a certain 
interesting set of rational functions. 

t Here Af(x, s) = f(x + 1, s) - f(x, s). 
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of a paper by Mrs. A. J. Pell appeared in the Bulletin of the American 
Mathematical Society, Ser. 2, Vol. 26 (1920), p. 149, indicating the character 
of this expansion problem for the case of two parameters. Consequently 
we shall leave out the formula? for this case. 

If we think of the several types of expansion problems — those for 
differential, difference, and integral equations — in intimate connection with 
the basic algebraic theory (developed in the memoir already referred to), 
it becomes apparent at once that the case of n parameters (for n > 1) is not 
confined to a set of n equations of the same sort. There is nothing to 
prevent one subset of the basic algebraic equations proceeding to differential 
equations as limiting forms, another to difference equations, another to 
integral equations, and so on. Thus we can see beforehand that we may 
formulate the expansion problem for a variety of mixed systems. In fact, 
we were led incidentally in § 2 to a special case of such mixed problems. 

Let us now indicate more precisely the nature of these mixed problems 
by considering the equations 

J n 

(!) -r 1 =H {an + X6# + ndj + vd i} )yj, i= 1, 2, • • •, n, 

ax j-i 

m 

(2) Asi(x) = ]T ifa + Xp<; + Wn + vTii)Sj{x), i= 1, 2, • • • , m, 

i=i 

(3) u{x) = \f a h K(x, $)u(f )df + pji L(x, £)«($)# + 

vf a b M{x, £)u{x, QdS, 

associating with them their adjoint equations 

dz- v-i 

(4) -j 1 = X) (- an - X6y f - ncji - vd H )zj, i= 1, 2, • • •, n, 
ax j=i 

m 

(5) At { (x) = J2 (— <Pn — ^Pn — Wh — vTji)tj{x +1), i = 1, 2, • • •, m, 

(6) ((*) = XX s K(i, x)v{£)dt + »f2 L(i, x)v(ij)di; + vf a b M(ij, *)»($)#. 
We suppose that the <p, p, a, r functions are analytic at infinity and vanish 
there to at least the second order. 

From (1) and (4), from (2) and (5), by methods already frequently 
illustrated, we have 

(7) Q/l2l + 2/2^2 + • • • + ynZn2xZa> = 0, 

m 

(8) £ [>,•( co )*,.( oo ) - Si(a)ti(a)2 = 0. 

We are now to set up homogeneous linear boundary conditions on y and on 
2 so that (7) is implied by them; and on s and on t so that (8) is implied 
by them. We shall suppose that this is done in such a way that the sets 
of characteristic values for the two problems [_(!), (2), (3) on the one hand 
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and (4), (5), (6) on the other] are the same and that the number of the sets 
is enumerably infinite. Then we shall denote the corresponding values, 
and solutions by 

X « M « „<*>, y w > z « ,.« ti w > „(*) > „w ; fe = 1, 2, 3, • • • . 

In equations (1), (2), (3) let us replace the symbols for the parameters 
and the functions sought by the corresponding symbols with the super- 
script k; in equations (4), (5), (6) let us introduce similarly the distinct 
superscript I; and let us denote in order the six new equations so formed 
by (1), (2), • • •, (6). From (1), (4) and the boundary conditions implying 
(7) we have readily 

n n n n 

(x<*> - x<o)y? EEW^+ (M w -M a) )// SSwW<fe 

(9) *= 1J=1 . .*= 1 '= 1 

+ („(*) _ v O> )f > £ £ d H y«W l) dx = 0. 
Similarly, from (2), (5) and the boundary conditions implying (8), we have 

oo m m 

(X<*> - X (i) ) ZEE P/<(« + w)*i(<* + w)tj(a + 1 + w) + ■■• 

/1AA W=0 i=l 7=1 

+ (*<» - i/») £ £ £ Tyi(« + «>)*(« + w)*X« + 1 + w) = 0. 

From (3) and (6) we get likewise 

nn (X<*> - \W)S a b Sa b Kti, x)tl<»(x)v«>Wdxdl; + • ■ ■ 

ui; + (*<*> - ^«)X 4 7? # (£, x)#(x) s ( ! '(9« = 0. 

Since the parametric differences X w — X (i \ • • •, involved in the last 
three equations are not simultaneously zero when h #= I the determinant 
of the coefficients of these differences must have the value zero when k 4= I- 
The relations thus arising express the generalized biorthogonality property 
of the solutions pertaining to the problem in hand. Once these are obtained 
the formal determination of the coefficients in the formal expansions 

00 

i = 1, 2, •••, n; j = 1, 2, •••, m, 

of a set of ran functions of three variables is easily effected by the method 
already employed several times in this paper. 

The foregoing we have outlined briefly as an instance of a sort of thing 
which may be done in a great variety of ways. The formal aspects of the 
problem are easily developed whether we have a single equation of a given 
type, as in the problem here treated, or any number of equations of each 
type. Also, we may employ in the list of types integro-difference and 
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integro-differential equations. In view of the developments of the section 
following it will be seen that partial differential and partial difference equa- 
tions are also admissible. Similarly one may allow equations involving 
functions of two variables and the operation of differentiation with respect 
to one and that of differencing with respect to the other; and so on. In this 
vast totality most types are without special interest; but a few actually 
arise in problems which come naturally to our attention. 

Nor is this all. It is also possible to treat similarly equations in which 
functions are subjected to different operations with respect to the same 
variable. We illustrate this with one of the simplest cases. Let /?* denote 
a set of n operations, the first k of them being identical with D and the 
remaining ones with A, where k is some one of the numbers 0, 1, • • •, n, 
and where D is the symbol for differentiation and A is the symbol for the 
operation Af(x) = f(x + 1) — /(x). Then form the system of equations 



(13) DiUi{x) = ^2 (fflfy + ^h)uj(x), i= 1,2, 



n, 



3= 



and associate with them the adjoint system 

(14) D{0i{x) = X (— a a — y>ji)vj(x + €j), i = 1, 2, •••,», 

where e» is zero or unity according as D t is D or A. If we multiply (1) by 
Vi(x + €i) and (2) by Ui(x), add the resulting equations member by member, 
and sum as to i from 1 to n, we have a result which may be put in the form 

n 

(15) ED i {u i (xMx)} = 0. 

Then, if appropriate conditions at infinity are satisfied, we have 

n 

^2 {m,-(=o )»,■(») — Ui{a)Vi{a)} = 0. 

It is clear now how one may complete the formulation of the expansion 
problem; the work proceeds in the closest analogy with special cases of it 
in §§ 2 and 4. 

6. Problems Arising from Partial Differential Equations and Integro- 
differential Equations. — As typical of the general situation in regard to 
expansion problems arising from the theory of partial differential equations, 
let us consider that which is to be associated with the equation 

(1) L{u) + \I*(u) = 

and its adjoint 

M{v) + \Mtiv) = 0, 
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where 



d 2 u , „, d 2 u , d 2 u , du . R du 
dx 2 dxdy dy 2 dx dy 



t / \ odu , du . 
Lxiu) = £ — + q— + ru, 
ox ay 

M («) = 7Tt ( fflB ) + 2 A" (M + Tl ( CT ) - T- M - I" W + ™ 
ox'' ctad?/ ajr ox dy 

M i(«) = — — - (pv) — — (qv) + ru, 
ox oy 

where the symbols u, v, a, b, c, a, /3, y, p, q, r denote functions of the in- 
dependent variables x and y. 

Let us multiply both members of equation (1) by v and of equation (2) 
by — u and add the resulting equations member by member; thus we have 

d f / du dv\ , , ( du dv\ , I . . da db\ 1 

rx[ a VVx- U d-x) + h Vd-y- U d-y) + { a + ^-dlc-ry) m \ 

. d V f du dv\ . , ( du dv\ 

+('+*-$rs)-]-°- 

Let S be any closed region in the xy-pl&ne whose boundary is cut twice 
and only twice by a line through it and parallel to either axis. Then if 
we integrate over S both members of equation (3) we have a relation 
which reduces readily to the form 

CV ( du dv\ . , / du dv\ , / , . da db\ "]"=* , 
... , r[ ( du dv\.i(du dv\ 

+(^ +x »-^-S)-]">-°' 

where r a and 5 a are limiting values of x in S for a given y and r and s are 
the extreme values of x and where p x , a x , p, a have like meanings under 
interchange of r61es of x and y. 

If in the first and second integrals in (4) we introduce the variable t 
in place of the given variable of integration by means of the relations 

t -V — P * _ x ~ r 
t = > t = > 

a — p s — r 

respectively, each integral reduces to an integral with respect to t from 
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to 1, and hence (4) itself may be put in the form 

(5) fo x Bdt = 

where B is a sort of bilinear form in the functions u and v and their first 
derivatives, the arguments of the functions being specialized in the way 
implied by the transformations involved (so that each of them is a function 
of a single variable t) . 

We now require linear homogeneous boundary conditions on u and v 
separately, independent of X, such that equation (5) is implied by the 
boundary conditions alone. We shall presently make clear the nature of 
these conditions by the explicit treatment of a particular case. For the 
present we shall suppose that these conditions have been set up in such 
way that the w-problem and the ©-problem have the same characteristic 
values X and that these are enumerably infinite in number. This done, we 
denote the characteristic values X and the corresponding solutions by Xi, 

Ui, Vif 1 =1 1 y Zy O, * * * . 

From (1) and (2) we now have the relations 

L(Ui) + XjXi(mj) = 0, 

M(v,) + \iMfr,) = (Xi - XjJJfifa). 

If we multiply the first of these equations member by member by v, and 
the second by — Ui, add the resulting equations member by member, and 
then integrate over 8, the resulting first member has the value zero on 
account of the boundary conditions, and hence the second member also' 
has the value zero. Omitting the non-zero factor X,- — Xy, this gives us 
the relation 

(6) ff Ui MM)dxdy = 0, i #= j. 

Similarly, we have 

(7) ffviLi{uj)dxdy = 0, i #= j. 

These two (equivalent) sets of relations express the generalized biorthogo- 
nality conditions for the problem in consideration. 

If we suppose that the integral f J*UiMi(vi)dxdy taken over 8 is dif- 
ferent from zero for every i we have a ready means of determining formally 
the coefficients ct in the formal expansion 

(8) f(x, y) = Y,c k u k (x, y) 

of a function f(x, y) of two variables in terms of the functions wj. In fact, 
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we obtain readily the relations 

fff(x, y)Mi(vk)dxdy 
(9) c k = ^ , &=1, 2, •••. 

ffu k {x, y)M 1 (v k )dxdy 

s 

In order both to make the method of procedure more explicit and to 
illustrate the greater simplicity of the self-adjoint case of the foregoing 
problem let us consider an instance associated with a class of equations 
of great importance in theoretical physics, namely, equations of the foi m 



d 2 u , d 2 u 
dx 2 dy 2 



(10) £ij+™+ X (0-M)«=.O, 



where g is a function of x alone and h is a function of y alone. For the 
region S we take now the square = x si 1, = z/ = 1; and we use the 
boundary conditions 

(n) u(t, 1) - u(t, 0) = 0, u(l, t) - u(0, t) = 0, 

u y (t, 1) - u y (t, 0) = 0, u x (l, t) - u x (0, t) = 0, 

where the subscript y or x denotes the partial derivative with respect to 
this variable. It is easy to verify that problem (10), (11) is self-adjoint; 
that is to say, if one writes the same problem in v one finds that these two 
problems, one in v and one in u, satisfy those conditions already imposed 
in the more general problems associated with (1) and (2). 
Equation (6) now has the special form 

(.12) fJ-fvig + h)UiVjdxdy =0, i 4= j, 

so that for the coefficients c k in (8) we have the simpler values 

(13) c k = fyfo 1 ^ + h ) u ^ x > y)f( x ' y) dxd v , k = 1 2 3 •■■ 

Jo 1 Jo 1 (g+ h){u k (x, y)} 2 dxdy 

Instead of considering all the solutions of (10) which satisfy conditions 
(11) we may fix our attention on any convenient class of them, as for in- 
stance those solutions which can be written in the form 

(14) uix, y) = P (x)<r{y), 

a class of particular solutions of great importance in theoretical and applied 
mechanics. Substituting u(x, y) of the form (14) into (10) and (11) we 
find that a constant n exists such that 

(15) p"(x) + \g(x)p(x) + ix P (x) = 0, p(l) - p(0) = 0, p'(l) - p'(0) = 0, 

(16) <T"(y)+M(y)<r(y) - ^(y) = 0, <r(l) - <r(0) = 0, <r'(l) - <r'(0) = 0. 
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Each of the last two differential systems is self-adjoint. Together they 
form [if written in a single variable x~] a system of two differential equa- 
tions (with boundary conditions) in two parameters X and p. General 
problems of this type have already been treated in § 3. It was through 
the emergence of the multiple parameter problem in connection with partial 
differential equations that I was first led to seek the formulation of the 
general multiple parameter problem for algebraic systems and also for 
the transcendental cases treated in this memoir. 

Let us derive the equations for the characteristic values of X and p in 
equations (15) and (16). Let pi(x, X, p) and p 2 (x, X, p) be a fundamental 
system of solutions of the differential equation in (15) and seek a solution 
p(x, X, p) so that the boundary conditions shall be satisfied. Constants 
Ci and c 2 are to exist such that 

p(x, X, p) = Cipi{x, X, p) + c 2 Pi(x, X, p); 

hence we have through use of the boundary conditions in (15) the relations 

Ci{pi (1, X, p) — Pi (0, X, ju)} + c 2 {p 2 (1, X, aO — Pi (0, X, At)} = 0, 
Cl { Pl '(l, X, m) ~ pi'(0, X, aO} + c 2 {p 2 '(l, X, aO - P*'(0, X, p)} = 0. 



The condition of consistency of these equations is that 



(17) 



Pi (L X, m) — Pi (0, X, At) P2 (1, X, /x) — Pi (0, X, aO 
pi'(l, X, aO ~ Pi'(0, X, aO P»'(l, X, it) - p 2 '(0, X, p) 



= 0, 



a condition on X and p alone. Similarly, if <ri(y, X, p) and <r 2 {y, X, p) are 
a fundamental system of solutions of the differential equation in (16), 
we are led to the following second condition on X and p : 



(18) 



<j\ (1, X, p) — o-i (0, X, p) <r 2 (1, X, p) — o- 2 (0, X, p) 
<n'(l, X, p) - 0-/(0, X, p) <r 2 '(l, X, p) - <r 2 '(0, X, p) 



= 0. 



The common solutions of equations (17) and (18) are the characteristic 
values for problem (15), (16). They are also the characteristic values for 
problem (10), (11) when the solutions u are restricted to be of the form given 
in (14). 

Let us consider for a moment the special case in which g = h = 1. 
Then equations (15) and (16) have solutions (among others) of the form 



where a and /3 are constants. The boundary conditions then require that 
we must have 

e a - 1 = 0, a(er - 1) = 0, e s - 1 = 0, ^ - 1) = 0. 
Hence a = 2miri and (3 = 2niri where m and n are integers, positive or 
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negative or zero. Then corresponding characteristic values X and m L"P er " 
haps not the complete set] are given by the relations X + ju = — 4mV 2 and 
X — m = — inV, or 

X = - 27r 2 (m 2 + n 2 ), ix = - 2?r 2 (m 2 - n 2 ), 

where m and n are integers. The corresponding solutions u(x, y) may now 
be written 

u mn (x, y) = e 2 " 1 *™^"™. 

The expansion of f(x, y) may be written 

+» 
f(x, V) = Y, c mn e 2mwix e 2nwi «. 



m, n=— oo 



(19) d3? 



For a real- valued function f(x, y) it is easy to transform this into an expan- 
sion in sines and cosines of integral multiples of x and y with real coefficients. 
Let us consider the extension of our problem to adjoint systems of partial 
differential equations involving a single parameter. As a typical instance 
we take the system 

d 2 Ui . d 2 Ui -Ak ( duj . , dui . \ 

+ df = £A aii dx- +K ^ +CijUi ) 

+ XX) [oiij -~ + Pa -^- + ytjuA , i = 1, 2, • • • , n, 

and its adjoint system 

dhi , dHi _ ^ [_ dvj , dvj ( _ da H _ db H \ 1 
(2Q) dx 2 + dy 2 ~ h V H dx bji dy + T V ~dx -dy-) Vj j 

Let us multiply the first of these by v, and the second by — «,-, add the 
resulting equations member by member, and sum as to i from 1 to n; thus 
we have 



(21) 



V l~— ( — * — — \ 4- —( ^ Ui — ^vA! 
tri[dx\ l dx ' dx) dy\ % dy * dy ) \ 

— Jl £ \q- I a a + ^«»7 [ •*¥>< + t- j hi + Xj8<y | •«,-«< = 0. 



After the method employed in the first part of this section we may now 
obtain from (21) the condition which must be realized as a consequence of 
linear homogeneous boundary conditions on the Ui and the t>,- separately. 
In (21) we integrate over a two-dimensional region .S of the xy-pl&ne, per- 
forming the integration with respect to x for those terms affected by the 
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operator d/dx and with respect to y for those affected by the operator d/dy. 
In this way we have a relation analogous to (4) above. Just as in the 
preceding case we introduce a variable t so as to reduce the two integrals 
to a single integral of the form 

(22) JVBdt = 

where B is a sort of bilinear form. The boundary conditions must be set 
up so as to realize this relation; they may be set up so as to realize the more 
special relation B = 0. 

Without attempting to characterize these boundary conditions in general 
[suggestions for forming them are afforded by' the earlier part of this 
section and by the treatment of a like problem in § 1] we suppose that 
they have been chosen so that those on the Ui are linear and homogeneous 
in the w* and their first derivatives and so that those on the v, are linear and 
homogeneous in the «,• and their first derivatives and that this has been 
done in such way as to realize condition (22) as a consequence of the bound- 
ary conditions. We shall suppose also that the characteristic values Xi, 
X2, X 3 , • • • of the M-problem are the same as those for the u-problem and 
that these values are infinite in number. The solutions corresponding to 
\k we denote by ui m , Vi W . 

In (19) let us write X fc and Ui W for X and u,; in (20) let us write X; and 
»/° for X and m, where k 4= I. Multiplying the first of these resulting 
equations by v^ and the second by — uf®, adding member by member, 
summing as to % from 1 to n, integrating as to x and y over the region S, 
simplifying by use of the relation to which (22) reduces when u t is replaced 
by «* w and m by i>,- (0 [a relation which is valid in view of boundary conditions 
of the sort described], we have a relation which becomes the following 
when the non-zero factor X/b — X; is dropped. 

(23) L* =1 i= l \ ox dy 

+ ( yii - d ^-^i) V ^ ) \ UiW ] dxdy=0 ' k * L 
Similarly we have the equivalent relation 

(24) //EE en, ^— + fa H£_ + y ifU w Vi v> dxdy =0, k * I. 

s L*=i i=i I dx dy J J 

Now let/f(a:, y), i = 1, 2, • • •, n, be a system of n functions of the two 
variables x and y and let us seek simultaneous expansions of the form 

(25) fi(x, y) = E cM"\x, y), i = 1, 2, • • -, n, 
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where the c„ are independent of i as well as of x and y. We seek to deter- 
mine the c„ on the supposition that the expression in the first member of 
(24) is different from zero whenever k = I. We multiply both members 

of (25) by 

du w du- w 

sum both members as to j from 1 to n and as to i from 1 to n, and integrate 
in the resulting equation as to x and y over the region S. In view of (24) 
it is seen that c* is the only coefficient left in the second member and that 
we have for C& the value 

■A v^ f du m du Xk) 1 

ff 2-, 2-, <*ij —£ h Pij -^— + 7ijUj W \ vt m (x, y)dxdy 

s i=i j=i [ dx dy j 

k= 1,2, ■■■. 

In case «*_,• = = fty relations (24) and (26) reduce to the notably 
simpler forms 



n n 



ffE E Jiju^hPdxdy = 0, k + l, 

S 1=1 j=\ 
n n 

ff Z) Z) VijUj W fidxdy 

c k = s i= : i= : k = 1,2,3, ■•-. 

ff Z) Z 7i 1 u 1 m '0i m dxdy 

S i=\ j=l 

Just as the problems in ordinary differential equations with one param- 
eter can be generalized, as we have seen, to the case of any finite number of 
parameters so the problems treated in the foregoing part of this section 
can be extended to partial differential equations and systems in any finite 
number of parameters. Moreover, any finite number of independent 
variables may be employed in place of the two variables x and y of the fore- 
going treatment. 

If we use the Volterra limiting process for passing from a finite number 
of differential equations to an integro-differential equation we are led from 
equations (19) and (20) above to certain adjoint partial integro-differential 
equations. For the sake of simplicity we shall confine our attention to the 
case in which the functions a„, &;_,-, an, /3;y are all identically zero. Then the 
integro-differential equations are the following : 

(27) dht(x,j, s) + d*u(x,j, s) = f/ Lc(x> y> s> f) + XTfe y> s> t)lu(x> y> ^ 
(2g) dh{x^y, s) + dh(x, y, s) = f ^ c{x> y> t> s) + Xy(x> y> t> s) y (x> y> t)dt 
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These equations may be treated by processes closely analogous to those 
by which (19) and (20) were treated. We summarize the main results. 
Corresponding to (23) and (24) we now have 

(29) ff[f* f/y(x, y, t, *)*<*> (x, y, <)«« (x, y, s)dsdt]dxdy= 0, k*l. 

s 

Then if we seek expansions of the form 

(30) f(x, y, s) = T,c v ^"\x,y,s) 
we have for the coefficients c k the values 

ffLf«f«y(x, V, s, t)u<- k) (x, y, t)f(x, y, s)dsdf\dxdy 

(31) Ck = SSLS/SJ'yix, y, s, <)«<*>(*, y, t)^{x, y, s)dsdQdxdy ' 

s 

k= 1,2, •••. 

7. On a Problem Arising in the Theory of Vibrating Plates. — An impor- 
tant differential equation in the theory of vibrating plates is the following: 

(1) d? +2 dx J df + dg i+ a di 2= °' 

Here a is a constant and £ is a function of x, y, t. In the theory of vibrating 
plates one desires to have a solution %(x, y, t) of (1) which f or t = reduces 
to a given function f(x, y) so that we have £(x, y, 0) = f(x, y). Such solu- 
tions $-{x, y, t) have been found in a variety of ways in terms of particular 
solutions of the form 

(2) $(x, y, t) = u(x, y)e int 

where n is a constant. If we substitute this value of £ into (1) we have for 
u the equation 

,3) i( „ ) + x„.g+2^L + g + X„-0, 

where X = — n 2 a. We shall formulate the expansion problem as con- 
nected with equation (3) in the general form* suggested by the standard 
procedure which we follow under the guidance of the basic algebraic theory 
which directs our investigation. 

The equation adjoint to the equation L{u) + \u = is obviously the 
equation L(v) + \(v) = 0. We have the following readily verified basic 
identity: 

* A number of special cases of this problem have been treated (with important results) 
in connection with the usual theory of vibrating plates. 
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t / \ _ r / \ _ d T d 3 u _ d 3 v _ dv d 2 u , du d 2 v 
dx\ dx s dx 3 dxdx 2 dxdx 2 



d 3 u . . . d 3 v dvd 2 u 

dxdy 2 dxdy 2 dxdy 2 

v 



dxdy 2 } dy\_ dy 3 dy 3 



dvd 2 u , dud 2 v 
dydy 2 dy dy 2 



*l .. . Li ? .2 !}.. Ait" 




dH . , _ 9 Jt d 2 u . 
da?dy dy dx 2 

In this identity s denotes an arbitrary parameter. 

Let S be any closed region of the :n/-plane whose boundary is cut twice 
and only twice* by a line parallel to either axis and passing through the 
region. Then if u and v are two solutions of (3) we have 

= ff\jL{u) - uL(v)2 dxdy = ff { £ Bi+ j- B 2 J dxdy, 

where Bi and B 2 denote the first and second brackets expressions in the 
second member of (4). From this we have readily the relation 

(5) fp'LBJ^dy + Sr S [B2lr p :dx = 

where r v , s v , r, s, p x , <j x , p, cr have the same meanings as in connection with 
equation (4) of § 6. By the same transformation as is employed in § 6 
we can reduce relation (5) to the form 

(6) fo'Bdt = 

where B is a sort in bilinear form in the functions u and v and their deriva- 
tives up to the third, the arguments of the functions being specialized in the 
way implied by the transformations involved (so that each of the functions 
is a function of the single variable t). 

We now require (as usual) linear homogeneous boundary conditions on 
u and v separately such that equation (6) is implied by the boundary condi- 
tions alone. A particular class of cases of great importance is that in which 
the boundary conditions imply the stronger relation B = 0. Leaving off 
for the moment the actual formulation of these boundary conditions let 
us suppose that they have been set up in such a way that the problem for 
the equation L(u) + \u = and the w-conditions has the same character- 
istic values X as the problem for the equation L(v) + X« = and the v- 
conditions and that these characteristic values are enumerably infinite in 

*The reader will see what modifications of the treatment are sufficient to deal with 
the case when these lines cut the boundary any finite number of times. 
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number. We denote the characteristic values X and the corresponding 
solutions by X,, U{, »,-, i = 1, 2, 3, • • • . Since equation (3) is self-adjoint it 
is seen that u t = vt for every i provided that the ^conditions and the v- 
conditions are the same. In this case the entire M-problem is self-adjoint. 
Now we have the relations 

L(u { ) + \im = 0, L{vj) + \n>j = (Xi — \j)v h i + j. 

Let us multiply the first of these equations by vj and the second by — u it 
add the resulting equations member by member, and then integrate over S; 
the resulting first member has the value zero on account of the boundary 
conditions, and hence the second member also has the value zero. Omitting 
the non-zero factor X» — Xy, this gives us the relation 

(7) ffu{i}jdxdy = 0, i #= j. 

s 

For the self-adjoint case this reduces to 

(7i) ffuiUjdxdy = 0, i =(= j. 

8 

If we suppose the integral in the first member of (7) to be different from 
zero when j = i (as it certainly is in the self-adjoint case) we have a ready 
means of determining the coefficients c& in the formal expansion 

(8) f(x, y) = 2 c k u k (x, y) 

k=\ 

of a function f(x, y) of two variables in terms of the functions u h . In fact 
we obtain readily the relations > 

fffix, y)v k (x, y)dxdy 
° k ~ J[fu k (x,y)vk(x,y)dxdy' ~ ' ' 

8 

In the self-adjoint case we have the simpler expressions 

Sff(x, y)Uk{x, y)dxdy 
(9l) Ck= S ff{u k {x,y)}Hxdy ' h=l > 2 >--'- 

8 

Suppose that we desire a solution %(x, y, t) of (1) which f or t = reduces 
to fix, y) so that we have formally 

CO 

£ix, y, 0) = X) CkUkix, y), 
where the c k have the values just obtained. Such a solution is formally 
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the following: 

£(x, y,t) = £ c k u k (x, y)e^ Ktl ' tt 

This is the actual solution desired provided that the boundary conditions 
are such as to render the solution unique and provided that the requisite 
conditions of convergence are satisfied. We have thus completed the 
formulation of the expansion problem for the present case. 

It is desirable to observe how easily one can extend the formal theory 
associated with equations (3) to (9) inclusive to that depending on the 
equation L(u) + \gu = in place of equation (3), where g is a function of 
x and y. The boundary conditions will be set up exactly as before, the 
function g not intervening in this part of the discussion. In the paragraph 
containing (7) we shall have a slight modification of formulae and there 
will result the generalized biorthogonality conditions 

ffguiVjdxdy = 0, i =¥ j; 

s 

and these give rise to the relations 

ffg(x, y)f(x, y)v k (x, y)dxdy 

1 — 1 O Q 

Ck ffg{x, y)u k {x, y)v k (x, y)dxdy ' ' ' ' 

s 

for determining the coefficients c k in the formal expansion (8). 

It is clear that the formal expansion theory connected with the equation 
L(u) + \gu = can be readily carried over to the system 

Z(«co) = £ (a i3 -u^ + \gifuP), i = 1, 2, ■ ■ ■, n, 

and to the simultaneous expansion of n functions in terms of the solutions 
obtained. This in turn can be made to yield heuristically the formulation 
of a corresponding problem in partial integro-differential equations. These 
two problems in one parameter can thence be easily extended to problems 
in any finite number of parameters. The method of making these exten- 
sions is sufficiently obvious in view of the methods already employed in this 
paper in closely related problems. 

In order to make clear the nature of the boundary conditions associated 
with equation (3) and described above through the intervention of equation 
(6), let us set them up explicitly and generally for the case when S is a 
rectangle and relation (6) is to be satisfied by requiring that B shall be 
identically equal to zero. We define S as that rectangle in which a Si x Si b, 
c Si y Si d. 
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Let us now write more explicitly 

(10) B 1 (x, y), 



B 2 {x, y), 



respectively, for the forms in the first and second brackets expressions in 
relation (4). Then for B we now have the value 



(11) 



B = (d - c){B 1 (b, c+ (d- c)t) - Bi(o, c + (d - c)t)} 

+ (b - a){B 2 (a + (b - a)t, d) - B 2 (a + (b - a)t, c)}. 



Now let us think of Bi(cc, y) as a sort of bilinear form in the two sets of 
variables 



du d 2 u d 2 u d s u d 3 u 



dv dH d 2 v d 3 v d 3 v 
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(12) u ' dx ' dx* ' dy 2 ' dx 3 ' dxdy 2 and *' dx ' dx> ' bf ' dx 3 ' dxdy 2 ' 

Let us form the matrix \\ai/.\ of this bilinear form where the element a^ 
of the ith column and /th row is the coefficient of the product of the ith 
function in the first set by the jth function in the second set. Then we have 



ffl»7 = 



It is easy to show that this matrix is of rank 4; for the last row may be 
made to consist of elements zero by adding to it (s — 2) times the second 
last row, the fourth row may be made to consist of elements zero by adding 
to it (s — 2) times the third row, the last column may be made to consist 
of elements zero by subtracting from it s times the second last column, and 
the fourth column may be made to consist of elements zero by subtracting 
from it s times the third column, the matrix so obtained being the following 
(which obviously is of rank 4) : 





1 



- 1 



In Bi(cc, y) we have the six variables u and the six variables v listed in 
(12). In B 2 (x, y) we have similar six variables u and six variables v. The 
matrix of the bilinear form B 2 (x, y) is of rank 4, as one may readily verify 
by the method used for Bi{x, y). Hence we may look upon B, defined as 
in (11), as a bilinear function of 24 variables u and 24 variables v, and this 
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bilinear form is of rank 16. Hence B may be put into the form 
(13) B = £ Ui{u)VM, 



where for each i Ui(u) is a linear homogeneous function of the 24 variables 
u and Vi(v) is a linear homogeneous function of the 24 variables v. The 
set Ui(u), i = 1, 2, • • •, 16, is linearly independent, and so also is the set 
Vi(y>), i = 1, 2, •••, 16. Then the boundary conditions which we take 
for the M-problem are 

(14) Ui(u) = 0, i= 1, 2, •••, 8, 

while those for the «>-problem are 

(15) Vi(v) = 0, i= 9, 10, •••, 16. 

|~For some parts of the theory one may take any number k of conditions 

Ui(u) = not greater than 15, i = 1, 2, • • •, k, and then use for the v- 

conditions Vi(v) = 0, i = k + 1, • • • , 16. But the form of conditions given 

in (14) and (15) is needful for certain parts of the development.] If the 

conditions on v in (15) are the same as those on u in (14) the problem is 

self-adjoint and the maximum of elegance is secured. 

University of Illinois, 
October, 1920. 



